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We study the system of two quantum dots lying on the central plane of a planar semiconductor microcavity. By solving 
the full Maxwell problem, we demonstrate that the rate of resonant excitation transfer between the two dots decays as 
dr 1 ! 2 as a function of the distance d at long distance. This very long-range mechanism is due to the leaky and guided 
modes of the microcavity, which act as effective radiative transfer channels. At short distance, the d~ 3 dependence of 
the Forster mechanism, induced by the electrostatic dipole-dipole interaction, is recovered. 
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1 Introduction Semiconductor quantum dots (QDs) 
are promising candidates for a solid-state implementation 
of quantum information technology. To this purpose, both 
orbital and spin electronic degrees of freedom might be 
exploited. One of the major challenges along this route 
consists in achieving a controlled quantum state prepara- 
tion and transfer of quantum excitation between different 
subsystems. This latter task is particularly difficult if the 
subsystems are spatially separated. Very effective control 
schemes might be based on controlled coupling of QD elec- 
tronic degrees of freedom with photons. Locally, strong 
coupling between an electron-hole pair state in a QD and 
a resonant mode of an optical cavity was recently demon- 
strated J!], as the semiconductor counterpart of atom-cavity 
strong coupling in cavity quantum electrodynamics. Thanks 
to optical control and detection, strong coupling should al- 
low full control of the QD quantum state. For QDs lying 
separate in space, the coupling to photons is the most ef- 
fective way to achieve controlled excitation transfer, partic- 
ularly by taking advantage of resonant photonic structures 
lIUGE]. The same mechanism could be used to transfer the 
electron spin state J5) . 

In order to get better insight into these mechanisms, the 
fundamental nature of the photon-mediated transfer must 
be studied. At short range, electrostatic coupling between 
the QD dipole moments results in the Forster resonant en- 



ergy transfer (5|, with a probability amplitude decaying 
as the third power of the QD distance. Recently, we have 
demonstrated Q[8) that the radiative field produces a long- 
range mechanism dominating over the Forster transfer at 
long distance. The mechanism can be pictured as the emis- 
sion of a photon by a QD and reabsorption by a distant QD. 
It decays as the inverse of the QD distance, as expected 
from spherical wave propagation, and could be exploited 
to engineer quantum correlations ITOl. Experimental evi- 
dence of it was recently found (9) . 

In this work we extend the analysis of radiative exci- 
tation transfer to the case where QDs lie in a planar semi- 
conductor microcavity (MC) with distributed Bragg reflec- 
tors (DBRs). In this configuration, the transfer takes place 
through both the main resonant mode of the MC and guided 
modes originating from the DBR structure. We present a 
semiclassical formalism in which the QD optical suscepti- 
bility is coupled to the full Maxwell equations. We prove 
that the probability amplitude of the spontaneous excita- 
tion transfer decays as d~ x l 2 as a function of the distance 
between QDs d. Depending on the direction of the electron- 
hole dipole moment and on the QD-cavity detuning, this 
ultralong-range transfer can be enhanced or inhibited, thanks 
to the resonant character of the MC. This principle might 
be exploited to control the transfer mechanism, e.g. by means 
of laser-controlled ultrafast cavity tuning ifTTl . 
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2 Theory We consider the system of two QDs em- 
bedded in a planar MC with DBRs. We assume a A-MC 
with the QDs lying on its central plane, where the elec- 
tric field has maximum amplitude. We assume the z-axis 
as the MC vertical axis, and the x-axis as the axis joining 
the two QDs. The method follows directly the one adopted 
in Ref . J7) . In the limit of low QD excitation, the excita- 
tion transfer is governed by linear optics and is described 
by the Maxwell equation for the electric field £ coupled to 
the linear optical susceptibility of the QDs 

,2 
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where r is the 3-D position vector, the dielectric con- 
stant of the MC central layer, and x me 3x3 linear optical 
susceptibility tensor of the QD subsystem. By assuming 
the strongest confinement of electrons and holes along the 
z-direction of the QD, we can assume that the same se- 
lection rules for the interband transition as in a quantum 
well hold. Then, by restricting to the electron-heavy-hole 
transition in a crystal of cubic symmetry (e.g. InGaAs), the 
susceptibility tensor is expressed as 
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where /i cv is the Bloch part of the interband dipole ma- 
trix element and \P a (r e , r/j ) is the electron-hole-pair en- 
velope wave function in the a-th QD. Given the strongly 
resonant character of the radiative transfer mechanism, we 
have safely restricted our analysis to the ground e-h pair 
state of each QD. For the QD wave function, we neglect 
Coulomb correlation and assume the factorized form: 



^ a (r e ,T h ) = fe(Ze)fh{zh)(/)e, a (Pe)(f>h, a (Ph) 



(3) 



For the calculations, the in-plane wave functions are mod- 
eled as isotropic Gauss functions with standard deviations 
a e and 07,. We further assume <5-confinement along the z- 
direction by taking fj(z) = S(z). This very simplified 
model can catch the essential features of the QD optical 
transition, relevant to the excitation transfer mechanism. 

Under these assumptions, the steps leading to a cou- 
pled mode equation are formally the same as in Ref. Q. 
Maxwell equation (Q]) ca be cast into an integral Dyson 
equation IT2l which finally reads 

E a ( W ) = K + £ GaP } u) E P ( U ) , (4) 

LU/3 — uo — lQ+ 

where Twj a is the transition energy of the a-th QD and 

E a («) = 5^a,kE*M , (5) 



(6) 



Here, ip a> k is the two-dimensional Fourier transform of 
ipa(p) = <t>e,a{p)4>h,a{p), which is assumed to be cen- 
tered at the position R a of the a-th QD on the plane, and 
Ek(w, z) is the Fourier-transform, with respect to the in- 
plane coordinate p, of the in-plane projection E(w,r) of 
the electric field £(oj, r). E° is defined as E Q and accounts 
for a possible external field E°. We have further defined 
ko = (oj / c)yfe~^. Eqs. (0]i-(|6]l are evaluated at z = 0, as 
the (5-assumption for the z-confinement allows to carry out 
the z-integral in (HJ explicitly. 

The central quantity of the present analysis is the in- 
plane Green's tensor of the photon Gk(co, z), which is de- 
fined as the Green's tensor of the two coupled equations 
for Ek(o;, z), resulting from Eq. (Q]i after having elimi- 
nated the uncoupled z-component of the electric field. No- 
tice that this quantity depends on k = |k| only, because of 
the circular invariance of the Maxwell problem. For QDs 
in a bulk dielectric medium, Gk(u), z) is given by a sim- 
ple analytical expression Q. In a MC, we have shown lTT3l 
that an analytical expression still holds. The tensor is diag- 
onal on the basis of transverse (TE) and longitudinal (TM) 
electric field polarization with respect to the direction of k. 
In particular, at z = we have 



GfcM = 2kl 
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where k z = \Jk^ — k 2 and r J k {uj) (j = t, b) are the com- 
plex reflection coefficients of the top and bottom DBR re- 
spectively, that can be computed by means of a transfer- 
matrix technique lfT4l . We define g^(^) and g£ (w) as the 
diagonal elements of Gfc(o>). They represent the propaga- 
tors for TM and TE modes respectively. In the limit r 3 k = 0, 
the in-plane propagator in a homogeneous medium Q is 
correctly recovered from Eq. ©. By replacing the (normal- 
ized) Gauss wave functions of the QD states and (O into 
©, straightforward algebra leads to a simple expression 
for the rates G a p(w). In particular, by assuming that the 
two QDs lie along the x-axis, the interdot tensor Gn (u) 
can be expressed in the (x, y)-basis as (for simplicity we 
omit the cj-dependence from the notation) 



M 



kd 



(9) 
(10) 



where d is the distance between the two QDs and a 1 = 
<T e <T ^/( <T e + a h)- ^1- ® indicates that, in the case of two 
isotropic QDs, the x- and y-components of the e-h polar- 
ization vector are separately conserved. This is in general 
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not true for an anisotropic QD, and the off-diagonal trans- 
fer rates will depend on the degree of anisotropy of the 
wave functions. 
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Figure 1 (a): Logarithmic grayscale plot of jf/J (w)|. (b): 
Same for \g%(u))\. 

3 Results and discussion We study the case of a 
GaAs MC with top (air) and bottom (substrate) DBRs made 
of 20 and 19.5 pairs of GaAs/AlAs A/4-layers respectively. 
We take n\ = 3 and n-i — 3.5 as the AlAs and GaAs re- 
fraction indices, and L c = 260 nm as the cavity layer thick- 
ness. We assume two identical QDs with a e = tr/j = 5 nm 
and square dipole moment p? cv = 0.66 meV • nm 3 , corre- 
sponding to a Kane energy of 21.5 eV as in In As. 

By inspection of Eqs. (IFll-dTOli, we can characterize the 
long-range behaviour of the radiative transfer mechanism. 
We first point out that Bessel functions J n (x) decay asymp- 
totically as x~ x l 2 . Hence, the J-2 contributions in (0 and 
( [Tol l dominate at long distance. As expected, if the dipole 
moments of the interband transition are aligned along the 
axis joining the two QDs (x -polarization), then the long- 
range transfer is mainly governed by longitudinal TM-modes, 
whereas for dipole moments orthogonal to this axis (y- 
polarization), it is the transverse TE-modes that dominate 
in the transfer mechanism. The Gauss function in $Q ex- 
tends to momenta much larger than fcrj, as a <C L c . Hence, 
computing the full mode structure of a DBR-MC is impor- 
tant in order to have a quantitatively good account of the 
transfer mechanism. In Fig.Q~]we display the absolute value 
of g^(iL>) and g^(uj)- Both quantities are characterized by 
one main cavity mode with a cutoff energy frujo ~ 1.36 
eV. At larger k, modulations due to the leaky modes of the 
DBRs are clearly visible. A main difference appear at large 
momenta, where the T-modes of Fig. \V[a) display another 
strong resonance corresponding to a guided mode close to 
fco, which is instead absent for L-modes. For k > Jcq, the 
DBR reflectivity rapidly approaches unity and the Green's 
functions approach those of a homogeneous medium. 

We turn to study the frequency dependence of the trans- 
fer amplitudes [G12 = x, y), which are plotted in 
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Figure 2 Transfer amplitudes |[Gi2(w)] X2; | (blue) and 
|[Gi2(w)] TO | (red) vs energy hid, at d = 300 nm. 

Fig.|2]for a distance d = 300 nm. They are strongly reso- 
nant with the main cavity mode at k = 0, indicating that 
the MC resonance still dominates in the radiative trans- 
fer process. The pronounced asymmetry, with significantly 
larger values at positive detuning, is due to the onset of 
the radiative transfer through the main cavity mode. Due 
to the transverse guided mode, [Gi2\ yy is about twice as 
large as [Gi2]i; 2 - when off resonance. At negative detun- 
ing the transfer is inhibited, although only by one decade 
because, as seen in Fig. Q] the contribution of leaky and 
guided modes is practically frequency-independent. This is 
very important in the light of proposals to use MC polari- 
tons as virtual (off-shell) intermediate states for spin-spin 
transfer 0. We show that resonant (on-shell) transfer may 
equally take place, independently of QD-MC detuning. 

The transfer amplitudes vs d are plotted in Fig. [3] for 
three values of Tllo. As a general trend, they show an os- 
cillatory behaviour, as the transfer occurs mainly via a few 
well-defined momenta given either by the main MC mode 
or by the guided mode. At very short distance all ampli- 
tudes display a d~ 3 -dependence. Our full Maxwell formal- 
ism thus correctly reproduces the Forster transfer mech- 
anism, caused by the electrostatic dipole-dipole interac- 
tion J6). Below the cavity resonance Hujq, the amplitude 
for x-polarized QDs decays as d~ 3 l 2 whereas the one for 
y-polarized QDs displays an ultralong spatial range char- 
acterized by a d _1 / 2 -decay. According to Fig. Q] in fact, 
the main contribution for x'-polarization comes from the 
transverse guided mode in g f k , which in Eq. © has an 
asymptotic dT^I 2 law. In Eq. ([Tol l instead, g\ multiplies the 
dominant long-range contribution Ji(kd) which decays as 
^-1/2 rp^g j ar g est ra diative transfer effect is obtained for 
u> = uiq, namely at resonance with the cavity mode. In 
this case, the resonant cavity mode is effective for both 
polarizations and we find a common d~ x l 2 decay. Quan- 
titatively, the rate is about 1 ns _1 at d = 300 nm, of the 
same order as the QD spontaneous radiative rate. The same 
qualitative trend is displayed when u) > uio, but the over- 
all transfer amplitudes are smaller, as expected also from 
the results of Fig. |2] Physically, the d~ x / 2 dependence is 
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Figure 3 Transfer amplitudes |[Gi2(w)]a; X | (blue) and 
|[Gi2(w)] TO | (red) vs QD distance d for (a): hu> = hu>o — 
0.03 eV; (b): hu = huj ; (c): hto = hui + 0.03 eV. The 
dashed lines are guides to the eye for the g?~ 3 / 2 dependence 
in (a) and for the d~ 3 and d" 1 ' 2 dependence in the three 
panels (see text). 



a direct result of cylindrical wave propagation in 2-D, and 
should be compared to the corresponding d^ 1 dependence 
for a bulk medium J7) . 

As in the case of QDs in a homogeneous medium [7, 
[U, the effect of radiative transfer is to produce collective 
many-QD eigenmodes - the polariton modes of the sys- 
tem. Their complex energies are the poles of the homoge- 
neous part of Eq. [4] (setting E° = 0). We plot these eigen- 
energies in Fig.|4]for the same three values of the QD-MC 
detuning as in Fig. [3] as a function of the (much smaller) 
detuning between the two QDs. As a general trend, the 
imaginary parts of the eigen energies (corresponding to the 
radiative rates of the collective modes) take different val- 
ues at zero QD detuning. One rate is much slower than for 
the uncoupled QDs, two take intermediate values and one 
is much faster. This collective effect is enhanced at reso- 
nance with the MC, as seen in the comparison of Figs. [4] 
(a)-(c). A similar conclusion is drawn for the real parts. 
These radiative effects are expected to become significant 
for a dense assembly of QDs embedded in a planar MC, 
either randomly or regularly distributed IT5l . or in the case 
of few QDs and light confinement in 3-D ffplUPl . 



4 Conclusions We have presented a semiclassical the- 
ory of radiative excitation transfer between QDs in a planar 
MC. Due to the 2-D nature of the system, the transfer am- 
plitude decays as d~ x / 2 with the QD distance. The trans- 
fer takes place through both the main cavity mode and the 
guided mode structure of the DBRs, with a rate compara- 
ble to the spontaneous emission rate. Additional light con- 
finement, e.g. through photonic crystal waveguides, would 
enhance the transfer rate, opening the way to applications 
in quantum information. 
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